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Abstract. We present a small data base of homogeneously derived photospheric element abundances of DZ white 
dwarfs and related objects. Our previous investigations are supplemented with the analysis of ultraviolet spectra 
for nine white dwarfs. Of particular interest is the detection of La absorption in van Maanen 2 and a determination 
of the effective temperature of this star. The new value is about 1000 K lower than previous results due to the 
strong ultraviolet absorption by metals which has to be considered consistently. The metal abundances of our 
sampl e stars are compatible with the predictions from the two-phase accretion model of Dupuis et al 



1992 



1993a| , [L993b| ). Small deviations can be observed for the abundance ratios in some objects. This could indicate 
non-solar metal-to-metal ratios in the accreted material. Hydrogen can be detected in virtually all of our objects. 
Ifowever, its average accretion rate must be at least two orders of magnitude lower than the metal accretion rate. 
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1. Introduction 

At the cool end of the white dwarf sequence traces of heav- 
ier elements can be found in several objects with hydrogen- 
and helium-rich atmospheres. The presence of metals is 
usually attributed to accretion from denser parts of the 
interstellar medium combined with diffusion into deeper 
layers. This idea has been put on a sound theoretical ba- 
sis by Dupuis et al. (1992, 1993a, 1993b| ) whose two-phase 
accretion model assumes that accretion of metals is low 
(M « 5 • lO"'^" M0 yr~^) most of the time and is larger 
(M w 5 • 10^^^ Mq yr^^) only during rare and short pas- 
sages through interstellar clouds. Since the accreted ele- 
ments diffuse downwards on time scales of typically 10^ 
years they are only visible during or shortly after such 
encounters. 
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Empirical tests of the accretion model benefit from dif- 
ferent diffusion time scales of the elements so that photo- 
spheric abundances and abundance ratios vary in a char- 
acteristic way with time after the encounter. Dupuis et 
al. (1993b) concluded that - within the large observa- 
tional uncertainties - the two-phase accretion model can 
account for the observed abundances and their ratios. The 
most notable remaining problem for the helium-rich ob- 
jects (spectral types DZ and DBZ) is the accretion of hy- 
drogen. As the lightest element, hydrogen is not subject to 
downward diffusion and should be accumulated in the up- 
per atmospheric layers. However, hydrogen is only rarely 
observed in DZ stars and the measured abundances are 
usually rather low. 

A general problem of previous comparisons is the un- 
certainty arising from observations with low signal-to- 
noise ratios and model atmospheres not including recent 
advances in input physics. In Paper I (Koester & Wolff 
pooo| ), we have concluded that it is worthwhile to study 
the abundance patterns in DZ white dwarfs again using 
new or newly calibrated observations and the latest im- 
provements in model atmosphere calculations. The ulti- 
mate goal is a homogeneous data base of element abun- 
dances for the comparison with the two-phase accretion 
model. 
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Table 1. Observations 



WU IN umber 


Name 


Observations 


A/ A 


AA/ 


0002+729 


GD 408 


lUE, SWP19006L 


1150-1950 


6.0 






lUE, LWR15058L 


1850-3250 


6.0 






MMT 


3750-4650 


1.0 


0038-226 


LHS1126 


HST/FOS, G190H 


1573-2330 


1.5 






HST/FOS, G270H 


2221-3301 


2.1 






HST/FOS, G400H 


3240-4822 


3.1 






BVRIJHK photometry 


4400-22000 




0046+051 


vMa2 


lUE, LWR03322L 


1850-3250 


6.0 






lUE, LWR06474L 


1850-3250 


6.0 






Calar Alto, 3.5m 


3730-5075 


4.0 






BVRIJHK photometry 


4400-22000 




0552-041 


LP 658-2 


HST/FOS, G190H 


1573-2330 


1.5 






HST/FOS, G270H 


2221-3301 


2.1 






HST/FOS, G400H 


3240-4822 


3.1 






BVRIJHK photometry 


4400-22000 




0752-676 


BPM 4729 


HST/FOS, G190H 


1573-2330 


1.5 






HST/FOS, G270H 


2221-3301 


2.1 






HST/FOS, G400H 


3240-4822 


3.1 






BVRIJHK photometry 


4400-22000 




1011+570 


GD 303 


lUE, SWP18994L 


1150-1950 


6.0 






lUE, LWR15047L 


1850-3250 


6.0 






MMT 


3750-4650 


1.0 


1225-079 


K 789-37 


lUE, SWP30001L 


1150-1950 


6.0 






lUE, LWP09836L 


1850-3250 


6.0 






SAAO, 1.9m 


3600-5000 


5.0 


1822+410 


GD 378 


lUE, SWP19549L 


1150-1950 


6.0 






lUE, LWR15578L 


1850-3250 


6.0 






MMT 


3750-4650 


1.0 


2216-657 


L 119-34 


lUE, SWP29720L 


1150-1950 


6.0 






lUE, LWR14743L 


1850-3250 


6.0 



Remarks 
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Bergeron et al. (1997) 



Bergeron et al. (1997 



Bergeron et al. (1997) 



Sion et al. (1988) 



Koester et al. (|1990|) 



Sion et al. (1988) 



We started in Paper I with the Hubble Space Telescope 
(HST) observations of the DZA white dwarfs L 745-46A 
and Ross 640. Both objects show weak Balmer lines in 
addition to metals and are therefore interesting with re- 
gard to the hydrogen problem mentioned above. We have 
used the high quality HST spectra to test our analysis 
methods and to study the remaining uncertainties of the 
model atmospheres. In comparison with previous inves- 
tigations, we were able to determine more accurate ele- 
ment abundances, especially for carbon and iron. Major 
improvements were possible for the modeling of the broad 
La lines in both stars which could be reproduced for the 
first time without any unphysical assumptions. This pro- 
vides new means for the analysis of the hydrogen content 
in cool white dwarfs. 

In this paper, we analyze all available ultraviolet obser- 
vations of white dwarfs with spectral types DZ and DBZ. 
The ultraviolet is very suited for the study of metal abun- 
dances because the relevant elements have in general much 
stronger lines in this part of the electromagnetic spectrum 
than in the optical part - with the exception of the Call 
H and K lines. We use for most stars observations with the 



International Ultraviolet Explorer (lUE) which have been 
newly calibrated (NEWSIPS) and which show significant 
improvements compared to older versions. HST/FOS ob- 
servations could also be used for three stars. 

2. Model atmospheres 

A general description of the procedures used for the calcu- 
lation of our LTE model atmospheres and synthetic spec- 
tra has been given by Finley et al. (1997). Further details 
can be found in Paper I. We only explain briefly the effects 
of different equations of states (EOS) and the treatment 
of ha absorption. 

We use for our models a version of the Hummer- 



Mihalas-Dappen EOS (Mihalas et al. 1990) which predicts 
larger non-ideal effects and larger ionization than the EOS 
by Saumon & Chabrier (1994) and Saumon et al. (1995). 
In Paper I, we have compared the results of our standard 
models with models were non-ideal effects are switched off 
and the ionization fraction is determined only by thermal 
effects using the unmodified Saha equation. For models 
with effective temperatures of about 8000 K non-ideal ef- 
fects are negligible but the situation may be different at 
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Fig. 1. Flux distribution of van Maanen 2: BVRIJHK photometry from Bergeron et al. (1997) and lUE observations 
(rebinned to 60 A resolution). The continuous curve is a model with T^ff — 5700 K, metals, and hydrogen, the dashed 
curve a pure helium model with Tcff = 6700 K 



the lower temperatures of some of the stars in this paper. 
We therefore test the influence of non-ideal effects in these 
cases. 

Of importance for the analysis of the HST spectra of 
L 745-46A and Ross 640 has been the La wing which dom- 
inates the UV spectra of both stars. Since model calcula- 
tions using van-der-Waals broadening by neutral helium 
result in a La wing extending far into the optical region 
we have developed in Paper I a new approach to model 
this line. Since we are only interested in the far wings 
we use the quasistatic limit of the broadening theory (see 
e.g. AUard & Kielkopf 1991) but take into account two 



effects which are important at high perturber densities: 
the variation of the transition probability with distance 
between emitter and perturber and the strongly repulsive 
potential at close distances. With this procedure we were 
able to reproduce the La wings in both stars at the same 
abundances as derived from optical Balmer lines. 

3. Observations 

Three of our objects have been observed with the Faint 
Object Spectrograph (FOS) onboard HST. For the re- 
maining six DZ/DBZ white dwarfs we have used low dis- 
persion lUE spectra. The lUE data were retrieved from 
the lUE final archive and have been calibrated with the 



NEWSIPS calibration (Nichols et al. |199J). In addition 
to these data, we have also used optical observations and 
optical photometric scans. A list of the observations can 
be found in Table 0. 

4. Analysis 

4.1. Van Maanen 2 

Van Maanen 2 (WD 00464-051) is one of the three clas- 
sic white dwarfs. It is a cool and very metal-rich DZ 



star, the optical spectrum exhibiting, besides the typi- 
cal Call H-hK lines, Cal at 4227 A, Fel at 4384 A, and 
several Fel lines on the blue wing of Call. Analysis of 
this star, using model atmospheres, has been pioneered 
by Weidemann (1960) and later by Wegner ( 1972 ) and 
Grenfell (1974). The ultraviolet spectrum is dominated 
by the strong Mg I line at 2852 A. Also important are the 
Mgll lines at 2796/2803 A and several weak Hues of Fel 
and Fe II. Until now, the ultraviolet spectrum has not been 
analyzed in detail. Bergeron et al. ( 1997[ [200 1| ) have re- 
cently determined T^s = 6770 K from BVRIJHK pho- 
tometry using a pure helium atmosphere. Combining this 
result with the trigonometric parallax they also derived 
log 5 = 8.40. 

We use the photometry together with the co-added 
lUE LWR spectra to define the flux distribution of 
van Maanen 2. An optical spectrum (3730-5075 A) ob- 
tained at the Calar Alto observatory complements the 
observations for the determination of metal abundances. 
We started the analysis with a pure helium atmosphere 
of Tcff = 6700 K. Fig. shows the flt of this model to 
the flux distribution. To account for the solid angle of 
the star the model has been scaled to reproduce the in- 
frared J magnitude. It fits the photometric data well but 
fails to reproduce the overall flux level at ultraviolet wave- 
lengths where the absorption from metals is important. 
Then we calculated model atmospheres and spectra with 
Teff = 6700 K and several chemical compositions to de- 
termine the abundances of magnesium, calcium, and iron 
from optical and ultraviolet lines by a visual comparison. 
It turned out that the backwarming effect from the ultra- 
violet absorption - especially from the strong Mg I line - 
and the effect of increased electron density significantly 
alter the flux distribution. The effective temperature has 
to be reduced and, consequently, the abundances have to 
be determined again. In an iterative process we finally de- 




Fig. 2. lUE spectrum of van Maanen 2 compared to a model with Tcff — 5700 K, \ogg = 7.9, logH/He = —5.0, 
logMg/He = —9.1, log Si/He = —9.0, logCa/He = —10.65, logFe/He = —9.4 (solid line) and a model with the same 
parameters but without hydrogen (dashed line) 



rived Tcff — 5700 K, about 1000 K lower than the original 
result from Bergeron et al. 

Fig. |l| shows that the new model can account for the 
complete observed flux distribution. At optical and in- 
frared wavelengths both model spectra are virtually iden- 
tical. Our effective te mperat ure is simi lar to previous re- 
sults fr om e. g. Wegner (|1972|) , GrenfeU (|l974[) , and Liebert 
et al. ( |l987|) who considered the contributions of metals 
as electron donators in their calculations. This shows the 
importance of an accurate chemical composition for de- 
termining effective temperatures of DZ white dwarfs and 
related objects - a result which is also found and discussed 



in more detail by Provencal et al. ( 2002 ) in their analysis 
of ProcyonB. 

From a visual comparison of the fits to the flux distri- 
bution with different models we estimate the error of Tcff 
to about ±200 K. The solid angle necessary to reproduce 
the photometry together with the trigonometric parallax 
from Hipparcos (tt = 226.95±5.35 mas; ESA |l997| ) and the 
mass-radius relation from Wood ( 1994 ; without hydrogen) 
implies a surface gravity of log f/ = 7.9. The uncertainties 
in Tefi and tt translate into ±0.2 for the log of the gravity. 
The higher val ue of logg — 8.4 derived by Bergeron et 
al. ( 1997 , 2001 ) is due to the higher effective temperature 
which gives a higher intrinsic flux and therefore a smaller 
radius and a larger mass for the star. 

The absorption from magnesium, calcium, and iron 
alone cannot reproduce the ultraviolet flux completely. As 
illustrated in Fig. ||, there is a significant discrepancy at 
A ^ 2300 A. We have included in the model calculations all 
lines from the identified metals and have also determined 
an upper limit for silicon. One may speculate that other 
unidentified metals may be responsible for the absorption. 
However, we do not think this to be plausible since, for ex- 
ample, the absorption from C I at 1930 A would result in a 
competely different shape and other metals than the ones 



considered here have not been observed so far in DZ white 
dwarfs. With regard to the experience with absorption 
from La in L 745-46A and Ross 640 (Paper I) we would 
rather attribute the missing opacity to hydrogen. Fig. |^ 
shows that adding an abundance of logH/He = — 5.0 im- 
proves the fit considerably. The Balmer lines are too weak 
at this abundance to be detected in existing optical spec- 
tra so that the analysis of the hydrogen content must rely 
entirely on La. 

The results for van Maanen 2 are summarized in 
Table ||. For the abundances, we include three different 
error sources: uncertainties from effective temperature, 
gravity, and from the choice of the equation of state. 
These systematic effects dominate the total error whereas 
pure statistical contributions are of minor importance. 
Formal statistical errors are comparitively low and do not 
reflect the true uncertainties. Therefore, we have deter- 
mined the individual contributions from the systematic er- 
ror sources using visual comparisons to models with differ- 
ent parameters and have added the contributions quadrat- 
ically. As mentioned in Sect. ^, non-ideal effects are not 
important in the somewhat hotter stars L 745-46A and 
Ross 640. However, a weak influence can be observed in 
van Maanen 2. If we use the Saha equation instead of our 
standard EOS, then the gas pressure increases whereas the 
electron density decreases and the absorption lines - espe- 
cially the wings - are in general somewhat broader. The 
abundances have to be changed by about 0.1 to 0.2 dex 
to account for this effect. 

4.2. LHS1126, LP 658-2, and BPM 4729 

The analysis of van Maanen 2 shows the importance of 
the ultraviolet region for the temperature determination 
and the question of the hydrogen abundance. It seems, 
therefore, worthwhile to investigate the ultraviolet spec- 
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tra of other white dwarfs with similar effective tempera- 
ture. In this section, we present HST/FOS observations 
of LHS1126, LP 658-2, and BPM4729, the coolest white 
dwarfs in our sample. The observations have been ob- 
tained in guaranteed time but have not been published 
so far since the fluxes differed decisively from the predic- 
tions of preexisting models raising some doubts about the 
accuracy of the FOS calibration. However, without go- 
ing into details, we consulted with Ronald Downes of the 
Space Telescope Science Institute and other members of 
the FOS team and have concluded that the data reduced 
with the standard FOS pipeline are ~ to the best of our 
knowledge - calibrated correctly. 

LHS 1126 (WD 0038-226) is one of the fe w C2 H stars, 
a spectral class introduced by Schmidt et al. ( 1995 ). These 
objects exhibit broad optical absorption features which 
have at first been interpreted as blueshifted C2 Swan 
bands. Schmidt et al., however, showed that C2II is the 
most probable carbon and hydrogen compound to be 
formed under the conditions in LHS 1126. The strong in- 
frared flux deficiency reported by Wickramasinghe et al. 
( |1982| ) could be explained by Bergeron et al. (1994) using 
collision-induced absorption of H2 with He. 

Bergeron et al. (1997, 2001) have determined Toff ~ 
5400±170K, logy = 7.91±0.17,andlogHe/H = 1.86 from 
optical/infrared photometry together with the trigono- 
metric parallax. In Fig. ^ we present a model calcula- 
tion (dashed line) using the parameters from Bergeron 
et al. and compare it with the HST/FOS spectrum and 
the photometry. The model flux is significantly too high 
at infrared wavelengths. This could be expected since we 
do not include collision-induced absorption in our calcu- 
lations. More severe is, however, the failure in the ultravi- 
olet. We have included absorption by La in the same way 
as for L 745-46 A and Ross 640 (see Sect. ^ and Paper I). 
The high hydrogen abundance of LHS 1126 results in an 
extremely strong absorption. 

The observed La profile can be fitted much better 
if the abundance is reduced to H/He — 3 • 10~^ (solid 
line in Fig. ^). La absorption is still significant as can be 
seen from the comparison with a pure helium atmosphere 
(dotted line). The low hydrogen abundance implied by 
the ultraviolet flux poses a severe problem since a high 
value is required to explain the infrared flux deficiency 
by collision-induced absorption. Our approach to model 
the La wing has been tested successfully at about 2000 K 
higher effective temperatures. It is possible that this ap- 
proach fails at lower temperatures although the observed 
profile can be fitted very well. It could nevertheless be 
expected that La absorption must be strong for an abun- 
dance of H/He « 10~^. However, the observations do not 
show a strong UV flux deficiency. 

A much simpler atmospheric composition is found in 
LP 658-2 (WD 0552-041) which is classified DZ due to th e 
presence of shar p Cal l lines (Eggen & Greenstein 1965|) . 
Bergeron et al. ( 2001 ) have fitted the optical and ultra- 
violet photometry using a pure helium atmosphere with 
Teff = 5060 ± 60 K and log 5 = 8.31 ± 0.02. In Fig. | we 



compare a model (dotted line) with these parameters to 
the observed flux distribution. The optical and infrared 
photometry can be reproduced but the ultraviolet part is 
again unsatisfactory. If the flux deficiency is attributed 
to the absorption from La then about H/He = 5 • 10"'' 
would be necessary (solid line in Fig. ^). In contrast with 
LHS 1126, the shape of the absorption feature cannot be 
reproduced in detail. Higher hydrogen abundances seem to 
be necessary for A < 3000 A but the absorption at longer 
wavelengths would then be too strong. Our approach for 
La seems to overestimate the absorption in the outermost 
wing. 

The third cool object with FOS observations is 
BPM 4729 (WD 0752-676). Wickramasinghe & Bessell 
( 1979|) have det ected weak Ha and H/3 lines. Bergeron et 
al. ( |l997' , ^001 ) could fit the optical/infrared photometry 
and Ha with Teff = 5730 ±110 K, log^ = 8.21 ±0.09, and 
a pure hydrogen atmosphere. We compare in Fig. || a pure 
hydrogen model (dotted line) to the observed flux distri- 
bution. Again, the ultraviolet part where La is important 
cannot be reproduced. The predicted proflle for La in a 
hydrogen atmosphere is completely at odds with the ob- 
servation. The ultraviolet part looks more like a helium- 
dominated atmosphere similar to LHS 1126 or LP 658-2. It 
is indeed possible to fit this region with H/He = 3 • 10^^ 
at about Teff = 5500 K. However, optical hydrogen lines 
are not visible under these conditions raising doubts on 
the validity of this solution. 

All three objects presented here show problems if the 
infrared and ultraviolet parts of the electromagnetic spec- 
trum are to be fitted simultaneously. In LP 658-2, addi- 
tional absorption, most probably from La, seems to be 
present which would not affect severely the solution from 
the optical/IR photometry. In LHS 1126 and BPM4729, 
however, the ultraviolet and infrared observations require 
completely different atmospheric compositions. We should 
note that the HST observations are consistent with exist- 
ing photometry in the B and U bands which is plotted as 
filled circles in Fig. |^. Therefore, the discrepancies between 
the observations and some models cannot be attributed to 
remaining uncertainties in the FOS calibration. 

Bergeron et al. ( 1997] , |200l| ) have also noted difficulties 
in fitting the blue energy distribution. They have tried to 
explain this with a missing opacity, attributed to a pseudo- 
continuum originating from the Lyman edge, for which no 
physical discription exists. However, one source of opacity 
which necessarily is present and which is very strong at 
cool temperatures is the La line. The ultraviolet part of 
the spectrum is therefore essential for a determination of 
the H/He ratio. We have shown with our new theoretical 
discription that at least in He-rich objects the wing ex- 
tends into the blue part of the optical spectrum. Under 
these extreme conditions our line profile calculations may 
be only a first approximation and we cannot yet present 
satisfactory fits. It seems to be important, however, to con- 
sider the whole energy distribution of these objects from 
the UV to the IR to achieve a correct interpretation. 
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Fig. 3. HST/FOS spectra (bars) of LHS 1126 (top), LP 658-2 (middle), and BPM4729 (bottom) com bined with 
VRIJHK photometry from Bergeron et al. ( |200l| ). Additi onal p hotometry in the B (Bergeron et al. 2001 ) and 
U band (averages from the compilation of McCook & Sion 1999| ) is plotted as filled circles. The observations are 
compared to varied model spectra. LHS 1126: T^s = 5400 K, \ogg = 7.9, H/He = 3 • 10~^ (solid line), pure helium 
(dotted), and H/He = 1 ■ 10"^ (dashed). LP 658-2: T^s 5050 K, log^ = 8.3, H/He = 5 • lO""* (solid), pure helium 
(dotted). BPM4729: Teg = 5500 K, logg = 8.21, H/He = 3 • 10"^ (solid), Tcff = 5730 K, logg = 8.21, pure hydrogen 
(dotted). All models have been fitted to the V magnitude 
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Fig. 4. La wing of L 119-34. The thick continuous curve is a model with T^ff = 9200 K and logH/He — —4.0. The 
dashed curve is a model with T^ff — 8900 K; the lower temperature can only partly be compensated by a lower 
hydrogen abundance (logH/He = —5.0). The observed spectrum shows also the CI line at 1930 A, CaH at 1840 A, 
and several Fe II lines 



Bergeron et al. ( 1997 , 2001) draw far-reaching conclu- 
sions based on their findings of H/He ratios in cool white 
dwarfs, e.g. the absence of He atmospheres between 5000 
and 6000 K, and their re- appearance below 5000 K. With 
our demonstration in this paper that the parameters from 
Bergeron et al. fail to explain the UV flux distribution 
of several objects it is clear that we have currently not 
reached a consistent interpretation of these objects. The 
question of H vs. He atmospheres or the H/He ratio in 
mixed atmospheres needs further study. 

4.3. L 119-34 and K 789-37 



L 119-34 (WD 2216-657) and K 789-37 (WD 1225-079) 
have similar effective temperatures of about 10 000 K. lUE 
spectra are available for both objects, for K 789-37 we 
can also use an optical spectrum (3600-5000 A) obtained 
with the 1.9m telescope of the South African Astronomical 
Observatory (Koester et al. 1990] ). K 789-37 exhibits Ha, 
H/3, Mgl, and Call in optical spectra. The lUE spectrum 
adds Mgll, Si I, Fel, and Fell. We use the lUE spectrum, 
the visual magnitude {V = 14.79) and the Stromgren col- 
ors {b-y = 0.05, u-b = 0.16) from Kilkenny ( |1986[ [1987D 
to define the overall flux distribution. 

We start the analysis of K 789-37 by investigating the 
effective temperature and the hydrogen abundance. Since 
La determines the spectral shape at the shortest UV 
wavelengths both quantities have to be analyzed simul- 
taneously. However, only a rough estimate of the hydro- 
gen abundances is necessary for the temperature deter- 
mination because La depends only weakly on the ex- 
act value of the abundance. If we assume logg = 8 we 
derive T^s — 10 500 K from a fit to the V magnitude 
and the lUE spectrum; a visual inspection of H/3 gives 
logH/He — —3.8. The Stromgren colors are also compat- 



ible with these values: b — y = 0.09 and u — & = 0.16 for 
Teff = 10 500 K. They would, however, favor a somewhat 
higher temperature of about 11 000 K. We can determine 
the surface gravity only very inaccurately with the ex- 
isting data. An estimate is possible due to the different 
sensitivity of H/3 and La to changes in log 5 and hydro- 
gen abundance. The effect of a lower or higher gravity on 
H/3 can be compensated to some extent by the hydrogen 
abundance but this leads eventually to an incompatible 
flux at the shortest lUE wavelengths. We estimate a pos- 
sible range of logg — 8.0 ± 0.5. The uncertainty of the 
gravity implies errors of ±300 K for Tcff and ±0.5 for 
logH/He. 

The new effective temperature is about 1000 K higher 
than the previous results from Liebert et al. ( |l987t 
Tes = 9700 K) using multichannel spectrophotometry and 
Koester et al. ( |199(]| ; T^e = 9500 ± 500 K) using lUE 
and optical spectra. The latter could reproduce the op- 
tical spectrum with 10 000 K whereas the UV flux of the 
model was too high with this temperature. We attribute 
these problems to the missing La opacity in their models. 

For the analysis of the metal lines, we use Tcs — 
10 500 K, logg = 8.0, and logH/He = -3.8. The de- 
rived abundances are listed in Tabic ^. The errors repre- 
sent the uncertainties in effective temperature and gravity. 
For magnesium, there are additional uncertainties due to 
the hue profile of the Mgll hues at 2796/2803 A. The 
observed profile cannot be reproduced exactly, a prob- 
lem which is also known from the analysis of Ross 640 
(Paper I) . The broadening of the magnesium lines is deter- 
mined by van-der-Waals interaction with neutral helium. 
To reproduce the observed profile, we have increased the 
broadening parameter for van-der-Waals inter action (75) 
by a factor of five (see also Zeidler-K.T. et al. 198(j ). We 
added (quadratically) an additional error of 0.5 dex to 
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account for the uncertainty due to the Une profile. The 
error of the iron abundance represents also the small dif- 
ferences if the abundances are derived from Fel or Fell, 
respectively. The silicon abundance has to be taken with 
caution because Si I cannot be identified clearly. The same 
is true for aluminum (log Al/He = — 7.75±0.5, not listed in 
Table From the non-detection of the C I line at 1931 A, 
we determine an upper limit of logC/He = —7.5. 

Contrary to K 789-37, hydrogen lines could not be de- 
tected in existing opti cal sp ectra of L 119-34, until now 
(see Zeidler-K.T. et al. 1986|) . However, the shorter wave- 
length part of the lUE spectrum indicates strongly the 
presence of La. We, therefore, determine simultaneously 
effective temperature and hydrogen abundance from the 
overall flux distribution as given by the lUE spectrum and 
the visual magnitude (V = 14.43, Eggen 1969). This is 
possible due to the different influence of T^s and La on the 
ultraviolet spectrum. The situation is illustrated in Fig.^. 
From a visual comparison we derive Tcs — 9200 ± 300 K 
and logH/He = — 4.015'q- Since the surface gravity can- 
not be determined with existing observations we have as- 
sumed log 5 = 8.0 ± 0.5. This is not critical because the 
shape of the ultraviolet spectrum depends only weakly on 
the gravity. The Stromgren colors b — y and u — b (Zcidler- 
K.T. et al. 1986 ) can also be used for a temperature es- 
timate. They give T^e = 9200^30°° K, in agreement with 
the above value. 



The element abundances, as determined from the lUE 
spectrum, are listed in Table ||. The errors reflect the un- 
certainties in effective temperature and gravity. For sili- 
con, an additional error of ±0.2 is added because of the 
difficulty to reproduce all Si I lines with the same abun- 
dance. An additional error of ±0.25 has also been taken 
into account for iron because of abundance differences be- 
tween Fel and Fell lines. 



4.4. GD408, GD378, and GD303 

GD408 (WD 0002+729), GD378 (WD 1822+410), and 
GD303 (WD 1011+570) are three DBZ stars with effec- 
tive temperatures well above 10 000 K. Photospheric Call 
lines could be detected in all three objects but other metal 
determinations do not exist in the literature. We examine 
t he lU E spectra and the optical observations of Sion et al. 

Weidemann & Koester ( |l99lD have determined Tes 
from multichannel spectrophotometry and \ogg with the 
help of the trigonometric parallax for CD 408. Oke et al. 
( |1984D used multichannel data to determine Tes and log g 
for CD 378 and CD 303. We have used the results for the 
gravity in our analysis but have tested the previous tem- 
perature determinations using the lUE spectra and the 
visual magnitudes. Our results for CD 408 and CD 303 
(see Table §) are compatible with the previous values; for 
CD 378 we deter mine a temperature about 1500 K hotter 
than Oke et al. (|l984|) . 



Mgll 

II 




2300 2400 2500 2600 2700 2800 
A / A 

Fig. 5. lUE spectra of CD 408 (top), CD 378 (middle), 
and CD 303 (bottom) compared with model atmospheres 
for the parameters of Table || 

The lUE spectra of all three stars reveal weak lines 
from Mgll and Fell (Fig. ||). Since these features are 
present in all three spectra we are confident that they are 
not observational artifacts. The derived abundances are 
listed in Table |2[ The errors include uncertainties due to 
errors in effective temperature and gravity. We ha ve also 
analyzed the optical spectra of Sion et al. (1988) again. 
Our new values for the calcium abundance are about two 
orders of magnitude higher than the values derived by 
Dupuis et al. (19931) because the scaling of the equiva- 
lent width by the square root of the abundance, as used 
by Dupuis et al. in their analysis, does not give reliable 
results here. 



Burkhard Wolff et al.: Element abundances in cool white dwarfs 



9 



Table 2. Effective temperatures, gravities, and element abundances (by number). The errors include systematic un- 
certainties and are determined from visual comparisons to model atmospheres. The five hottest objects are of spectral 
type DBZ(A) and the others of spectral type DZ(A) 



Object 




Tcff/K 


log 3 


logH/He 


logC/He 


log Mg/He 


log Si/He 


log Ca/He 


logFe/He 


vMa2'' 


1 


5700±200 


7.9±0.2 


-5.0±0.5 




-9.1±0.3 


< -9.0 


-10.65±0.2 


-9.4±0.2 


L 745-46A'' 


2 


7500±200 


8.0±0.2 


-3.1±0.4 


< -8.0 


-9.05±0.2 


-9.4±0.35 


-10.60±0.2 


-9.8±0.2 


Ross 640'' 


3 


8500±200 


8.0±0.2 


-3.3±0.3 


< -9.0 


-7.25±0.65 


-7.5±0.5 


-8.65±0.45 


-8.3±0.15 


L 119-34" 


4 


9200±300 


8.0±0.5 


-4.01?- 


-6.0±0.5 


-6.8±0.7 


-7.9±0.4 


-9.1^* 


-7.95±0.7 


K 789-37" 


5 


10500±300 


8.0±0.5 


-3.8±0.5 


< -7.5 


-7.6±0.6 


-7.5±0.5 


-7.9±0.2 


-7.4±0.3 


GD 408" 


6 


13750±250 


8.0±0.2'' 


-6.0±0.3'' 


< -7.0'' 


-8.5±0.5 


< -8.0 


-9.6±0.3 


7 c+0.5 
— ' -0-1.0 


HS2253" 


7 


14700 ±500 


8.0±0.25 


-5.0l«- 


< -7.5 


-5.4±0.2 


-5.5±0.7 


-6.1±0.4 


-5.1±0.1 


GD 40^= 


8 


15150±600 


8.0±0.25 


-5.3±0.5 


-7.2±0.2 


-6.0±0.6 


-6.5±0.5 


-6.7±0.3 


-6.2±0.1 


GD 378" 


9 


17000±800 


7.9±0.2-^ 


-4.0'' 


< -6.7'' 


-6.5±0.8 


< -6.5 


-8.15±0.5 


< -5.5 


GD 303" 


10 


18000±1000 


7.8±0.2^' 


< -5.5" 


< -6.5'' 


-5.75±0.5 


< -7.0 


-7.75±0.5 





this paper; '' Paper I (Koester & Wolff |2000| ) ; " Fr iedrich et al. ( |l999|); 
Weidemann & Koester (|l99l[); ^ Oke et al. (|l984|); » Sion et al. ([1988D 



* Zeidler-K.T. et al. 019861 ): 
'' Wegner & Nelan (p87|) 



5. Discussion 

We have analyzed effective temperatures, gravities, and 
surface abundances for six DZ/DBZ white dwarfs. For the 
following discussion, we add the results from Pape r I an d 
the two DBZ white dwarfs from Friedrich et al. ( |l999|) . 
These ten stars span virtually the complete range of tem- 
peratures of DZ/DBZ white dwarfs. They are the only ob- 
jects with ultraviolet observations which provide the best 
possibility to investigate the whole spectrum of metals; in 
optical observations, only Call is found in general. We 
collect the parameters for all stars in Table ||. 

We start with a comparison of the observed abun- 
dances with the pr edicted ranges from the model of 
Dupuis et al. ( 19931j ) in Fig. |[ The two continuous curves 
in each panel represent the lowest and highest possi- 
ble abundances at each effective temperature. The lower 
curves correspond to the equilibrium abundances reached 
in the low accretion phase (M « 5 • 10-^*^71/0 yr~^ for 
49 million years in the model of Dupuis et al.) whereas 
the upper curves correspond to the maximum abundance 
reached during passages through interstellar clouds (M « 
5 ■ 10""'^^ A/q yr"^ for one million years). Both curves 
should be taken with some caution since the actual values 
of the accretion rates are not known exactly. Moreover, 
the upper curve depends also partly on the duration of 
the cloud encounter since an equilibrium between accre- 
tion and diffusion is not reached at all temperatures during 
the assumed passage value of 10^ years. 

The comparison shows that the model assumptions 
can account in general for the observed metal abundances 
since most observed values are well between the two limits 
of the model predictions. Only the observations of GD 40 
and HS 2253-^8023 (numbers 7 and 8 in Fig. ||) are in 
part above the upper limits but this can be attributed to 
the mentioned uncertainties. The data in Fig. ^ seem also 
to indicate a correlation of the abundances with effective 



temperature. This is, however, a selection effect since at 
higher temperatures only high metal abundances can be 
detected in the low resolution and low signal-to-noise ob- 
servations normally used to identify DZ white dwarfs (see 
e.g. Zeidler-K.T. et al. |1986 for quantitative estimates of 
the visibility ranges). 

Further information about the accretion process can be 
extracted from the ratios of the metal abundances. Fig. |^ 
shows the values for the six independent combinations of 
the four observed metals. Again, the continuous curves 
define the allowed ranges from the model of Dupuis et 
al. ( 1993b ). The straight line is the solar abundance ratio 
which is reached in the first moments of a cloud encounter 
if the accreted material has solar composition. Afterwards, 
the ratio is altered due to diflFusion and approaches the 
equilibrium value (not shown in Fig. ^ which is always 
within 0.5 dex of the solar ratio. The second curve in Fig. |^ 
corresponds to the maximum (or minimum) value reached 
shortly after a cloud encounter when the abundance of the 
element with the shorter diffusion time scale is again near 
the value during low accretion. Afterwards, the ratio turns 
back to the equilibrium value. At any time, the observed 
abundance ratio must be between the two curves shown 
in Fig. 0. The difference between them depends on the 
diffusion time scales of the two elements. 

Most individual abundance ratios fall into the pre- 
dicted range. They are also often near the solar value 
(and the equlibrium value) which can be seen well in the 
Mg/Ca and Fe/Mg diagrams. This implies that we observe 
the DZ stars mostly during the high accretion phase - in 
agreement with the observation that the abundances in 
Fig. 1^ are near the upper curve. The deviations from the 
predicted regions could be attributed to individual varia- 
tions of the chemical composition of the accreted material 
which could be different from star to star. Howeve r, we 
follow here the discussion of Dupuis et al. (1993b) who 



explained these differences by the same changes of the 
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Fig. 6. Observed metal abundances (error bars) compared with the predicted ranges from the two-phase accretion 
model of Dupuis et al. ( 1993b| ). The lower continuous curve is the equilibrium abundance reached in the low accretion 
phase. The upper curve is the maximum abundance during an encounter with an interstellar cloud. Observational 
upper limits are denoted by dashed bars. The labels correspond to the numbering in Table |^ 



composition for all stars. This conclusion is supported by 
systematic tendencies which were not clearly present in 
previous studies: Several observed values especially of the 
Si/Ca and Si/Fe ratios are below the allowed range. 

If the abundance of one element in the accreted mate- 
rial is changed then the whole region of possible ratios in 
Fig. is shifted by that amount. If we want virtually all 
measured values to fall into the predicted range then the 
silicon abundance must be changed by a factor of «0.2 



(—0.7 dex) and the calcium abundance by a factor of ?a2 
(0.3 dex). Similar values were also determined by Dupuis 
et al. ( |l993b| ). 

The underabundance of silicon is interesting in connec- 
tion with the assumption that metals are accreted onto 
white dwarfs in the form of grains. Silicon has the lowest 
condensation temperature of t he fou r metals usually found 
in DZ spectra (see e.g. Spitzer 1978 , Fig. 9.1). During the 
formation of grains for example in the winds of cool gi- 
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Fig. 7. Observed metal abundance ratios (error bars) compared with the predicted ratios from Dupuis et al. (1993b). 
The two continuous curves give the possible range if the accreted material has solar composition. Upper limits from 
the observations are denoted by dashed bars 



ants it condenses later than the other elements. It has, 
therefore, a lower depletion factor in the gas phase of the 
interstellar medium and may be underabundant in the 
grains accreted onto white dwarfs. During the infall onto 
the white dwarf silicon could be evaporated earlier than 
other metals and may partly not reach the surface. Other 
elements with higher condensation temperatures would 
survive the accretion process better. Note that calcium, 
which is probably overabundant, has the highest conden- 
sation temperature of the four metals. 



In contrast to heavier elements, hydrogen is not ac- 
creted in the form of grains. As mentioned in the in- 
troduction, it can only be observed in DZ white dwarfs 
at unexpectedly low abundances. Nine of our objects ex- 
hibit hydrogen, in two of them we could show for the first 
time its presence through absorption by La. The distri- 
bution of hydrogen abundances against the effective tem- 
peratures of the stars is shown in Fig. ^. As in previ- 
ous investigations, a correlation with effective tempera- 
ture cannot be observed. We show as a comparison also 
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Fig. 8. Observed hydrogen abundances compared to the 
expected abundance for continuous accretion of material 
with solar composition with M = 10~^^ Mq yr^^. 



Fig. 9. Calcium-to-hydrogen ratios compared to the max- 
imum Ca/H ratio which could be expected for accretion 
of material with solar composition. 



the expected hydrogen abundance for continuous accre- 
tion of M = IQ-^^ Moyr-i (Dupuis et al. |l993bD . The 
calculation assumes solar composition of the accreted ma- 
terial and a start of the accretion at Tcff = 20 000 K. The 
exact choice of the starting point is not critical because 
of the rapidly increasing cooling ages towards lower tem- 
peratures. Fig. H shows that the actual average accretion 
rate of hydrogen must be lower than M = 10^^^ M© yr~^ 
for our sample objects. On the other hand, the accretion 
rate for metals ~ averaged over 5 • 10^ years - in the model 



our improved data set than in the previous investigation 



of Dupuis et al. is about M = 10 



-16 



M, 



©yr 



Therefore, 



the hydrogen accretion rate is at least two orders of mag- 
nitude lower than the metal accretion rate, even for those 
objects with the highest hydrogen abundances. 

In this context, the ratio of hydrogen to other metals 
is also interesting. Fig. ^ shows the observed calcium-to- 
hydrogen ratios. The weak tendency in temperature which 
might be observed in the data is due to a similar effect 
visible in the calcium abundances (see above). Fig. ^ shows 
also the maximum value of Ca/H which can be expected 
for accretion of material with solar composition within 
the model of Dupuis et al. ( 1993bD . All observed values 
are significantly higher than the theoretical upper limit. 
This shows again that hydrogen is not accreted at solar 
abundance ratios relative to metals. 



of Dupuis et al. (1993b). Hydrogen is present in virtu- 
ally all of our objects. However, the average accretion rate 
of hydrogen must be more than two orders of magnitude 
lower than the average accretion rate of the metals - even 
in those objects with the highest hydrogen abundances. 
Therefore, the fate of hydrogen remains unclear. 

Our success in modeling La absorption in L 745-46A 
and Ross 640 (Paper I) provides a new means of detect- 
ing hydrogen in those cases where the Balmer lines are al- 
ready very weak. This method has been used to determine 
the hydrogen abundance in L 119-34 and van Maanen 2. 
This is important especially for the latter object because 
Ha is not visible at this combination of effective tempera- 
ture and hydrogen abundance. In three other white dwarfs 
with temperatures similar to van Maanen 2 the ultravio- 
let region of the spectrum indicates a completely differ- 
ent chemical composition than that implied by the optical 
and infrared photometry. However, we are at the moment 
not in the position to solve this discrepancy. The analysis 
of van Maanen 2 shows also that the analysis of optical 
and infrared photometry can lead to wrong effective tem- 
peratures if the absorption by metals in the ultraviolet is 
not considered. The temperature would be about 1000 K 
higher without appropriate blanketing in the UV. 



6. Conclusions 

The element abundances of the ten DZ and DBZ white 
dwarfs presented in this paper form a homogeneous set of 
data which has been compared to the two-phase accretion 
model of Dupuis et al. (1993b). The absolute values are 
in the expected range for accretion from the interstellar 
medium. Most metal abundance ratios are also compat- 
ible with the predictions. Small deviations could be an 
indication of non-solar metal-to-metal ratios in the ac- 
creted material. This tendency is more clearly visible in 
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